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We present an analysis of bulk and surface electronic structures of thallium based ternary III-V- 
VI2 series of compounds TIMQ2, where M=Bi or Sb and Q=S, Se or Te, using the ab initio density 
functional theory framework. Based on parity analysis and (111) surface electronic structure, we 
predict TlSbSe2, TlSbTe2, TlBiSe2 and TlBiTe2 to be non-trivial topological insulators with a single 
Dirac cone at the T-point, and TlSbS2 and TlBiS2 to be trivial band insulators. Our predicted 
topological phases agree well with available angle-resolved photoemission spectroscopy (ARPES) 
measurements, in particular the topological phase changes between TlBiSe2 and TlBiS2- Moreover, 
we propose that Weyl semimetal can be realized at the topological critical point in TlBi(Si- x Se x )2 
and TlBi(Si- cc Te cc )2 alloys by breaking the inversion symmetry in the layer by layer growth in the 
order of Tl-Se(Te)-Bi-S, yielding six Dirac cones centered along the Y — L directions in the bulk 
band structure. 
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I. INTRODUCTION 

Topological insulators are a new class of materials, 
which have attracted intense interest in the last few 
yearsi"— due to their exotic properties. These materi- 
als support an odd number of surface state bands with 
linear dispersion in the bulk energy gap, which can be 
viewed as a sea of massless Dirac Fermions. The conduct- 
ing surface states in topological insulators are protected 
by time-reversal symmetry and are immune to scatter- 
ing by non-magnetic impurities, thus opening new av- 
enues for dissipationless transport. There also is a vig- 
orous ongoing search for topological superconductors^, 
with the possibility of realizing Majorana fermions, which 
are their own antiparticles 6 with potential application to 
quantum computing 7 . Topological insulators have also 
generated a considerable excitement due to the possibil- 
ity of exploring the Higgs mechanism and realization of 
a Weyl semimetal in a condensed matter systemA"— 

A topological phase was initially predicted in 2D 
HgTe/CdTe quantum wellsi^ and subsequently verified 
experimentally 15 . Soon thereafter a 2D conducting sur- 
face state was realized in the bulk band gap of three di- 
mensional thermoelectricsi^"— Bii-^Sb^, Bi2Se3, Bi2Te3 
and Sb2Te3. However, magnetotransport studies^ 7 
have shown that bulk transport dominates in these ma- 
terials, motivating continued search for other 3D topo- 
logical insulators with single Dirac cone surface states 
residing in the bulk energy gap. Since then first- 
principles calculations have suggested a large variety of 
topologically interesting materials ranging from oxides 22 
to the Heusler family of compounds^ 3 - - — . Another 
class is thallium based ternary III-V-VI2 chalcogenides, 
which were proposed theoretically 26 -^ 7 - and then verified 
experimental!}* 2 ^"— . These studies showed the existence 



of single Dirac cone type surface states at the T-point 
in TlSbSe 2 , TlSbTe 2 , TlBiSe 2 and TlBiTe 2 . Also, it has 
been found that p doped TlBiTe2 superconducts^i, where 
superconductivity is attributed to six leaf like bulk pock- 
ets in the Fermi surface and the surface state becomes 
superconducting 29 . 

Interestingly, recent studie a 32 i 33 of TlBi(Si_ x Se :c )2 al- 
loys show that a topological phase transition can be re- 
alized by modulating either spin-orbit coupling or the 
crystal structure, and that the surface Dirac fermion be- 
comes massive at the quantum phase transition. This 
can be viewed as a condensed matter version of the Higgs 
mechanism in which a particle acquires a mass by sponta- 
neous symmetry breaking. This system thus may provide 
a model system which connects condensed-matter physics 
to particle physics. It has also been proposed that a Weyl 
semimetal phas o 34 i 35 could be achieved at the phase tran- 
sition between a topological and normal insulator if we 
explicitly break time reversal symmetry 8,9 or inversion 
symmetr y 10 ! 11 . In this new phase, valence and conduc- 
tion bands touch at certain points, called Weyl points, 
where dispersion is linear. These Weyl points come in 
pairs with positive and negative helicities and are robust 
to perturbations in the bulk material. 

Our motivation for undertaking the present study is 
to provide a comprehensive investigation of the bulk 
and surface electronic structures of the thallium based 
ternary III-V-VI2 series of compounds TIMQ2, where 
M=Bi or Sb and Q=S, Se or Te, within a uniform 
first principles framework. In particular, not only the 
nontrivial compounds are studied 36,37 but the topolog- 
ically trivial compounds TlBi(Sb)S2 is included since 
the recent ARPES measurement observed the topolog- 
ical phase transition in the alloy TlBi(Se,S)2 between 
nontrivial TlBiSe2 and trivial TlBiS2- 33 We compare our 
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theoretical predictions with the available ARPES results. 
Insight into the topological nature of these compounds is 
gained through slab computations in which the thickness 
of the slab is varied. Finally, the possibility of realizing a 
Weyl semimetal phase through strain in TlBi(Si_ x Se ;E )2 
and TlBi(Si_ x Tea;)2 alloys is examined for the first time, 
and it is shown that, when the inversion symmetry of the 
system is broken, a Weyl phase is indeed possible at the 
topological critical point where the system undergoes a 
transition from a trivial to a non-trivial (topological) in- 
sulator. Bulk and surface electronic structures of many 
members of the six ordered Tl-compounds considered in 
this work have been discussed in several earlier papers in 
the literature^^^. 

The organization this article is as follow. Section II 
gives details of the bulk and surface computations. In 
section III, we explain the bulk crystal structures, elec- 
tronic structure and the parity analysis used to infer the 
topological nature of various compounds. Section IV dis- 
cusses the slab structure used, relaxation and size de- 
pendent effects. In section V, we explain the topological 
phase transition and realization of the Weyl semimetal 
phase in TlB^Si-^Se^ and TlBi(Si_ x Te :c )2 alloys. Sec- 
tion VI summarizes conclusions of this study. 



II. COMPUTATIONAL DETAILS 

Electronic structure calculations were carried out 
within the framework of the density functional theory 38 
using VASP 39 (Vienna ah initio Simulation Package), 
with projected augmented wave basis 40 . Exchange- 
correlation effects were treated using a generalized gradi- 
ent approximation^ and the spin-orbit coupling (SOC) 
effects are included as implemented in the VASP package. 
For structure optimization a plane wave cut-off energy of 
350 eV and a T-centered 8x8x8 k-mesh with conjugate 
gradient algorithm(CGA) 4 ^was used. Lattice parameters 
and ionic positions were adjusted until all components of 
Hellman-Feynman force on each ion were less than 0.001 
eV/A. Our surface electronic structure calculations are 
based on a slab geometry using bulk relaxed parameters 
with a net vacuum of 15 A, a plane wave cut-off energy 
of 350 eV, and a T-centered 9x9x1 k-mesh. Since ionic 
relaxation s 36 ! 37 are important, ionic positions in all slabs 
were optimized until the z-component of the Hellman- 
Feynman forces was less than 0.005 eV/A. Topological 
phase transition in TlB^Si-^Se^ and TlB^Si-^Te^ 
systems was investigated by taking the critical concen- 
tration of sulfur to be x=0.5 33 . 



III. RESULTS AND DISCUSSION 
A. Bulk Crystal Structure 

Thallium based ternary chalcogenides III-V-VI2 share 
a rhombohedral crystal structure (space group Dg d (R3 



m)) with four atoms per unit cell which occupy the Wyck- 
off positions 3a, 3b, and 6c 43 . We illustrate the crys- 
tal structure with the example of TlBiTe2, which can be 
viewed as a distorted NaCl structure with four atoms in 
the primitive unit cell 4 ^ - — and a sequence of hexagonal 
close packed layers in the order Tl-Te-Bi-Te (see Fig.l). 
The conventional unit cell is hexagonal and contains 12 




FIG. 1. Primitive (A) and conventional (B) hexagonal crystal 
structure for TlBiTe2 . The conventional (hexagonal) unit cell 
has three times more atoms than primitive cell.(C) ELF plot 
of TlBiTe2 (see text). (D) Brillouin zone for the primitive 
unit cell with four time reversal invariant points T, F, Z and 
L. 

atoms. The hexagonal lattice constants an and ch can 
be computed from the rhombohedral parameters via the 
relation 



an = 2a# sin — , ch — clr V3 + 6 cos a (1) 

where a# is the rhombohedral lattice constant and a 
is the rhombohedral angle. Notably, the present com- 
pounds form a layered structure similar to binary topo- 
logical insulators A 2 B 3 (A=Bi, Sb and B=Se, Te) 16 
where atoms are arranged in a sequence of five atomic 
layers or a quintuple repeating pattern. The bonding 
within the quintuple of layers is ionic- covalent type while 
between one quintuple group and the next is van-der 
Waals type 4 ^. In order to gain insight into the nature of 
bonding in the thallium based chalcogenides we have per- 
formed an ELF (Electron Localization Function) 47 study, 
which shows that there is a negligible value of the ELF 
between the thallium and tellurium layers, indicating pre- 
dominantly ionic type bonding between these layers. In 
contrast, the value of ELF between the bismuth and tel- 
lurium layers is significant implying ionic- covalent type 
bonding. The optimized lattice and internal parameters 
u for the six Tl-compounds considered are given in Table 
1. 



TABLE I. Optimized lattice constant clr, angle a and internal 
parameter u for six thallium based rhombohedral compounds. 



Compound 


CLR 


a 


u 




(A) 


(degrees) 




TlSbSe 2 


7.955 


30.869 


0.2374 


TlSbTe 2 


8.421 


31.367 


0.2386 


TlSbS 2 


7.718 


30.568 


0.2352 


TlBiSe 2 


7.989 


31.368 


0.2391 


TlBiTe 2 


8.367 


31.847 


0.2413 


TlBiS 2 


7.979 


29.764 


0.2368 



B. Bulk Band Structures 

The bulk band structures of TlSbQ2 Q=S, Se, and Te, 
shown in Fig 2, indicate these three compounds to be 
narrow gap semiconductors. TlSbSe2 is seen to be a di- 
rect band gap semiconductor with the valence band max- 
imum (VBM) and conduction band minimum (CBM) ly- 
ing along T — L. Here the spin-orbit coupling plays an 
important role as it induces a band inversion at T sug- 
gesting a possible nontrivial topological phase, a point 
to which we return below. TlSbTe2 and TlSbS2 are also 
seen from Figs. 2(B)-(C) to be direct band gap semicon- 
ductors with VBM and CBM at the T-point. Since Te is 
heavier than Se, there is a large spin-orbit (SOC) effect 
in TlSbTe2 . But S atom is lighter than Se and Te so that 
SOC effects are smaller in TlSbS2- 

The bulk band structures of TlBiQ2 (Q=Se, Te, and 
S) are shown in Fig. 3. There is a large effect of spin- 
orbit coupling due to presence of the heavier Bi atom. 
Our band structures agree well with an earlier study 45 , 
except in the case of TlBiTe2 , which was previously found 
to be 45 an indirect gap semiconductor with the VBM 
lying along the L — Z direction and the CBM at the 
T-point. However, for our relaxed structure (including 
spin-orbit coupling) the VBM and CVM still lie at the 
aforementioned ^-points, but the system is semi- metallic 
with a band gap of -10 meV. This is in accord with the 
corresponding experimental results 2 ^, which found the 
material to be semi-metallic with a band gap of -20 meV. 




C. Parity analysis 

Since all the investigated compounds possess inversion 
symmetry, a parity analysis 4 ^ can be used to identify 
the Z2 topological phases. There are eight time reversal 
invariant points in the rhombohedral Brillouin zone but 
only four points (T, F, L and Z; see Fig. 1 (D)) are 
inequivalent. Products of parity eigenvalues at these four 
momenta are given in Table II with and without spin- 
orbit coupling. 

Table II, shows that product of parity eigenvalues in 
TlBiSe 2 , TlBiTe 2 , TlSbSe 2 and TlSbTe 2 changes at T 
as spin-orbit coupling is turned on, yielding a nontrivial 
topological invariant Z2 = 1. The non-trivial topological 



FIG. 2. Bulk band structures of rhombohedral (A) TlSbSe 2 
(B) TlSbTe 2 (C) TlSbS 2 along high symmetry lines with (red 
lines) and without (green lines) spin-orbit coupling. 



character of these compounds is thus due to band inver- 
sion at the r-poin t 26 i 27 i 36 . On the other hand, there is 
no band-inversion for TlSbS2 and TlBiS2 at any of the 
time reversal invariant momenta, indicating that these 
compounds are topologically trivial. 



4 








-1 



-2- 



V 

■'J TIBi 




v V 








tt? (ft 












^ 



FIG. 3. Bulk band structures of rhombohedral (A) TlBiSe 2 
(B) TlBiTe 2 (C) TlBiS 2 along high symmetry line with (red 
lines) and without (green lines) spin-orbit coupling. 



IV. SURFACE ANALYSIS 



Compound 


r 


Xl 


L 


x3 


F 


x3 


z 


Xl 


z 2 




wso 


SO 


wso 


so 


wso 


so 


wso 


so 




TlSbSe 2 


+ 


— 


+ 


+ 


+ 


+ 


+ 


+ 


1 




_l_ 




_l_ 


_l_ 


_l_ 


_l_ 


_l_ 


_l_ 


I 


TlSbS 2 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 





TlBiSe 2 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


1 


TlBiTe 2 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


1 


TlBiS 2 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 






TABLE II. Products of parity eigenvalues at the four inequiv- 
alent time-reversal invariant k-points in the six investigated 
compounds, 'wso' refers to without spin-orbit coupling and 
'so' to with spin-orbit coupling. 



minations, depending on which atom lies in the topmost 
layer. The bond length between Tl-Te is large (d=350 
pm), whereas bond length between the Bi-Te is small 
(d=318 pm). Thus, out of the four possible surface ter- 
minations, we have used the one with Q (Se,Te and S) 
atom at the surface and Bi (Sb) under it in the second 
layer as this termination has a minimum number of dan- 
gling bonds 26 . Optimized bulk parameters were used to 
construct slabs of different thicknesses. A symmetric 39 
atomic layer slab of TlBiTe2 and the associated surface 
Brillouin zone is shown in Fig. 4. Surface relaxation is 
known to play an important role 36 ' 37 , and accordingly, 
we have relaxed all atomic layers in the slabs used, al- 
though the relaxation effect is small as we move into the 
bulk. 



B. TlSbQ 2 (Q=Se,Te and S) 

The Q-atom terminated surface electronic structures 
of the three compounds are shown in Fig. 5 for different 
slab thicknesses together with the associated projected 
bulk bands(blue colored region). Figs. 5(A) and (B) 
show the surface band structure of TlSbSe2 for slabs of 
35 (thickness ~ 6.4 nm) and 47 (thickness « 8.6 nm) 
layers, respectively. In the 35 layer slab there is a band 
gap of 50 meV at the F-point. As the number of layers 
increases, the size of the gap decreases. For 47 layers, this 
gap is negligible and we get a clear Dirac-cone surface 
state in the bulk gap region. Moreover, around -0.8 eV 
we obtain a Rasbha-type, trivial spin split surface state 
in both slabs 36 . Fig. 5(C) shows results for a 47 layer 
TlSbTe2 slab, which are similar to those for TlSbSe2- 
Finally, TlSbS2 does not display any surface state, which 
confirms its topologically trivial nature. 



Slab Structure and Relaxation 



C. TlBiQ 2 (Q=Se, Te and S) 



A hexagonal unit cell is used for all slab computations 
with atomic layers stacked in the z-direction. Since atoms 
along the z-direction are sequenced in the order Tl-Te (S, 
Se)-Bi(Sb)-Te(S, Se), there are four possible surface ter- 



Similar to TlSbQ2, TlBiQ2 also shows thickness depen- 
dent slab electronic structure. Fig. 6 shows the surface 
band structure of TlBiSe2- With a slab thickness of ^6.3 
nm we get a gapped surface state with a gap of 36 meV 
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FIG. 4. (A) Structure of a 39 atomic layer slab of TlBiTe 2 
stacked along the ^-direction. (B) 2D Brillouin zone with 
three special k-points F,M and K. 
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FIG. 5. Surface electronic structure of various slabs: (A) 35 
layer («6.4 nm) TlSbSe 2 ; (B) 47 layer (« 8.6 nm) TlSbSe 2 ; 
(C)47 layer (« 9.2 nm) TlSbTe 2 ; and (D)39 layer (« 6.9 nm) 
TlSbS2, along high symmetry lines in the surface Brillouin 
zone. The bulk bands projected on the surface Brillouin zone 
are shown in blue color. 



and as we increase slab thickness to ~7.0 nm this gap 
reduces to 9 meV. The gap at small slab thickness arises 
due to interaction between the opposite faces of the slab. 
As we further increase the slab thickness to « 8.7 nm 
and ^10.9 nm we obtain a clear Dirac state with a neg- 
ligible gap. At a critical thickness of ~ 7.0 nm the Dirac 
point is isolated, i.e. no other states are at the energy of 
the Dirac point, which agrees with recent angle-resolved 
photo-emission spectroscopy (ARPES)2S^^ measure- 
ments. The observed Dirac state for slabs thicker than 
7.0 nm is surrounded by the surface bands, which is some- 
what different from ARPES results and requires further 
study. 

Fig 7(A) shows results for TlBiTe2, which exhibits a 
Dirac-like surface state in the bulk gap region at the 
thickness of ~ 7.4 nm (39 layers). The Dirac point lies 
0.1 eV below the Fermi energy in accord with ARPES 
studies 29 . The electronic structure displays an indirect 
band gap near the f point with bulk conduction band 
minimum at the f -point and bulk valence band maxima 
along the T — M and f — K directions. The location 
of Dirac point in our calculation is close to experimental 
results, where Dirac point lies at 0.3 eV below the Fermi 
energy. Results for TlBiS2 in Fig. 7(B) do not have any 
metallic surface state in the bulk gap region indicating 
its topologically trivial character. 



V. TOPOLOGICAL PHASE TRANSITION AND 
WEYL SEMIMETAL 

A. TlBi(Si_ x Se x ) 2 

The topological phase transition was observed in 
TlBi(Si_ x Se x )2 in which the S and Se atoms are disor- 
dered and inversion symmetry is preserved on the aver- 
age with the inversion center at the (Se,S) site^ At the 
critical composition x~0.5, the conduction and valence 
bands meet at the Y point and form 3+ ID Dirac-cone 
bulk states. The Dirac-cone states are doubly degenerate 
due to inversion and time-reversal symmetries. When the 
inversion symmetry is broken, the spin degeneracy can be 
lifted. Here we consider an ordered phase of TIBiSSe with 
layers in the order Tl-Se-Bi-S as shown in Fig. 8(E) and 
(F). The inversion symmetry is now seen to be broken 
as Se and S sites are no longer inversion centers. Our 
theoretical fully relaxed structure, which is rhombohe- 
dral, is found to be stable, indicating that the material 
should be possible to realize via molecular beam epitaxy 
(MBE) techniques. In the absence of inversion symmetry, 
parity analysis cannot be applied to delineate the topo- 
logical nature, but adiabatic continuity arguments can 
be used. Accordingly, we start from the normal insula- 
tor with a large value of c/a and systematically reduce 
its value. During this process the system passes from 
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FIG. 6. Surface electronic structure of various slabs: (A) 39 
layer (« 7.0nm) TlBiSe 2 ; (B) 35 layer (« 6.3nm) TlBiSe 2 ; and 
(C)47 layer (« 8.7nm) TlBiSe2, along high symmetry lines in 
the surface Brillouin zone. Projected bulk bands are shown 
in blue color. 



normal to a topological insulator as shown in the Fig. 
8(A). The band structure for the normal and topological 
phase through the critical point is shown in Figs. 8(B)- 
(D). The critical point occurs when the gap closes. At 
c/a=4.98, the band gap is 0.135 eV, and the bulk valence 
band has no s-character while the bulk conduction band 
possesses a finite s-character at T. These characteristic 
symmetries are the same as in TlBiSe2, indicating that 
the two compounds are adiabatically connected and are 
thus both topologically nontrivial. With increasing c/a 
value the gap decreases and becomes zero at the critical 
value c/a=5.33. Upon further increasing c/a, the bulk 
gap increases and becomes 0.129 eV at c/a=5.58 and the 
bulk valence and conduction bands swap their orbital 
character at the T-point. The system is now adiabati- 
cally connected to TlBiS2 and is, therefore, topologically 
trivial. 

Interestingly, at the topological critical point, the band 
gap at T remains finite with a value of 0.062 eV. How- 
ever, the band gap closes along r — L at fc=(0.0, 0.0, 
0.007). Thus, by breaking inversion symmetry, we obtain 
a nondegenerate spin-polarized bulk Dirac cone, instead 
of the doubly degenerate Dirac cone found in the Se/S 
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FIG. 7. Surface electronic structure of: (A) 39 layer (« 
7.4nm) TlBiTe 2 slab, and (B) 47 layer (« 8.6nm) TlBiS 2 slab, 
along high symmetry lines in the surface Brillouin zone. The 
blue color shows projected bulk bands. 



disordered TlBi(Se,S)2- The bulk valence and conduction 
bands now touch each other and display linear dispersion, 
which can be described by a two component wave func- 
tions, i.e. by the Weyl equation. A Weyl semimetal is 
formed at the critical point with six Weyl points centered 
along r — The detailed band structure is shown in 
Fig. 9. 



B. TlBi(Si_.Te.) 2 

Analysis of TlBi(Si_ a; Te ;E )2 follows along the lines of 
that of TlBi(Si_ ;c Se a ;)2 in the preceding subsection, and 
therefore, we only make a few relevant remarks. Here 
again we take the critical point to be at x=0.5, and 
break the inversion symmetry by considering an ordered 
phase with layers in the sequence Tl-Te-Bi-S, and com- 
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FIG. 9. Band dispersion at the critical point in 
^ TlBi(Si-xSe x )2. Inset shows dispersion near the T-point. 



ting is larger than in TlB^Si-^Se^. 



Nontrivial 



FIG. 8. (A) Phase diagram for the topological phase tran- 
sition. The system is tuned from the normal to topological 
insulator by changing the c/a ratio. Panels (B)-(D) show the 
bulk band structure of TIBiSeS for c/a=5.22, c/a=5.33 and 
c/a=5.58, respectively. Primitive rhombohedral (E) and con- 
ventional hexagonal (F) crystal structure for TIBiSeS. The 
layered structure is with no inversion symmetry. 



pute the band structure for a series of c/a values. The 
results, summarized in Fig. 10, show that at c/a=4.93 
there is a direct band gap of 0.123 eV along F — L, 
although there is an indirect gap of ~ lOmeV. At 
the critical value of c/a=5.00, the gap becomes zero 
along r - L at k = (0.0,0.0,0.025), even though the 
gap remain finite at T with a value of 0.210 eV. With 
further increase in c/a, the gap reopens and becomes 
0.221 eV at c/a=5.57. Adiabatic continuity arguments 
combined with the change in the orbital character of 
the bulk valence and conduction bands at the critical 
value of c/a then allow us to conclude the non-trivial 
to trivial insulator transition shown in Fig. 10(A). As 
to the Weyl semimetal phase, at the critical c/a value, 
TlBi(Si_ x Te :c )2 displays six Weyl points along the Y — L 
direction at k = (0.0,0.0,0.025), although the spin split- 
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FIG. 10. (A) Phase diagram of topological phase transition 
in TIBiTeS. Panels (B)-(D) show the bulk band structure of 
TIBiTeS for c/a=4.93, c/a=5.00 and c/a=5.57, respectively 



VI. CONCLUSIONS 

We have carried out an ab initio study of bulk and sur- 
face electronic structures of six thallium based III-V-VI2 
ternary chalcogenides TIMQ2, where M(Bi, Sb) and Q(S, 
Se or Te), with focus on delineating the topological na- 
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ture of these compounds. TlBiTe2 is found to be a semi- 
metal with a band gap of -10 meV while the other five 
Tl-compounds are all small band gap semiconductors. 
Based on an analysis of parities of bulk band structures, 
we predict that TlSbSe 2 , TlSbTe 2 , TlBiSe 2 and TlBiTe 2 
are non-trivial topological insulators with band inversion 
at the T-point, but TlSbS 2 and TlBiS 2 are trivial band 
insulators. Moreover, surface state computations show 
that the surface Dirac states lie in the gap region at the T- 
point in all four aforementioned topological compounds. 
Our predicted topological phases and the T-point cen- 
tered Dirac-cone surface states are in substantial accord 
with available ARPES results. Electronic structures of 
slabs with different numbers of layers were computed in 
order to gain insight into thickness-dependent effects. 
The gap opens at the Dirac point for thin slabs and 
decreases with increasing thickness. Finally, we inves- 
tigated TlBi(Si_ x Se :c ) 2 and TlBi(Si_ :c Te :c ) 2 alloys for 
£=0.5 where the inversion symmetry was explicitly bro- 
ken by using layers in the sequence Tl-Se(Te)-Bi-S and 
the c/a ratio was varied. Both alloys were found to un- 
dergo a topological transition at a critical value of c/a 



at which the spin degeneracy of the Dirac states is lifted 
and a Weyl semimetal phase could be realized with six 
Weyl points in the bulk Brillouin zone located along the 
r — L directions. 
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